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Abstract 

We examine a spatial evolution of lightwaves in a nonlin- 
ear photonic crystal with a quadratic nonlinearity when 
simultaneously a second harmonic and a sum-frequency 
generation are quasi-phase-matched. We find the condi- 
tions of a transition to Hamiltonian chaos for different 
amplitudes of lightwaves at the boundary of the crystal. 
PACS: 42.65.Sf,05.45.Ac,42.70.Mp,42.65.Ky 

Wave mixing in nonlinear optical materials is a basis 
of modern optical sciences and technologies. Cascad- 
ing several wave-mixing processes in the same low-loss 
material one can in principle achieve a high efficiency 
using a large value of the lowest-order optical nonlin- 
earity. The theoretical investigations of cascading of 
several scalar optical three- wave-mixing processes in the 
bulk materials with x^ nonlinearity has a long history 
0. In particular, Akhmanov and co-workers have found 
the efficiency of a third harmonic generation (THG) via 
cascading of a second harmonic generation (SHG) and 
a sum-frequency mixing (SFM) in a quadratic medium 
[Q, while Komissarova and Sukhorukov have described 
an efficient parametric amplification at a high-frequency 
pump in the same system || . Obviously, the observation 
of these nonlinear effects demands a simultaneous sat- 
isfaction of phase-matching conditions for several para- 
metric processes as perfectly as possible. On other hand, 
it has been shown later that the systems, for which sev- 
eral optical wave-mixing processes can be simultaneously 
phase-matched, are in general nonintegrable; therefore a 
competition of two (or more) parametric processes can 
often result in a chaotic spatial evolution of lightwaves 
[ij. However, until nowadays it was unclear how to 
achieve a phase-matching for several processes in homo- 
geneous medium employing traditional techniques, such 
as using a birefringence in ferroelectric crystals. 

The solution of this problem has been found rather re- 
cently 0, D ; it consists in an introduction of the differ- 
ent types of artificial periodicity of a nonlinear medium 
resulting in a formation of nonlinear ID and 2D super- 
structures termed optical superlattices B or nonlinear 
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photonic crystals (NPCs) In NPCs there is a pe- 

riodic (or quasiperiodic) spatial variation of a nonlinear 
susceptibility tensor while a linear susceptibility tensor 
is constant. 

In these engineering nonlinear materials a phase mis- 
match between the interacting lightwaves could be com- 
pensated by the Bragg vector of NPC. The idea of such 
kind of quasi-phase-matching (QPM) was introduced by 
Bloembergen and co-workers many years ago fl~l|j . How- 
ever, only recently the rapid progress in a fabrication 
of high quality ferroelectric crystals with a periodic do- 
main inversion made the QPM method very popular 
[|[ [l2]| . We should stress that the conditions for QPM 
may be fulfilled for several wave-mixing processes simul- 
taneously; the QPM also has an advantage as using of 
largest nonlinear coefficient. 

Nowadays there are several experiments on an obser- 
vation of third and fourth harmonics in different peri- 
odically or quasiperiodically poled ferroelectric crystals 
with nonlinearity J?], [y| H, which clearly demon- 
strate an importance of multiple-mixing in NPCs for the 
potential applications. Modern theoretical activities on 
the nonlinear lightwaves interactions in NPCs are mainly 
focused on the studies of strong energy interchange be- 
tween the waves Jl2| (this is a development of the earlier 
activities @ |), as well as on the formation of spatial 
optical solitons |ll| . 

In this work we describe a novel for the physics of 
NPCs effect of Hamiltonian optical chaos. Namely, we 
show that spatial evolution of three light waves, par- 
ticipating simultaneously in SHG and SFM in the con- 
ditions of QPM, is chaotic for many values of complex 
amplitude of the waves at the boundary of x (2) -NPC. 
There also exists an integrable limit when the evolution 
of waves is always regular regardless of absolute values 
of their complex amplitudes. The integrable limit corre- 
sponds to the particular values of two combinations of 
wave phases at the boundary of nonlinear medium. In 
particular, the problem of THG belongs to the integrable 
limit, therefore in the conditions of recent experiments 
10, [D| nonlinear light dynamics should be always 
regular. However, even a rather small change in ampli- 
tudes and phases of waves at the boundary of crystal, in 
respect to those considered in |7], [T^, [l4|], should result 
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in a transition to chaos. 

We consider a spatial evolution of three co- 
propagating plane waves 
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E = — Aj exp [i(jujt — kj 
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where index n takes only odd values. Now we substi- 
tute this expansion into Eqs. ([!]), take into account the 
QPM conditions (|^) and make an averaging of resulting 
motion equations over "the short characteristic spatial 
scale" 27r/A. We have the following basic equations 



in a periodically poled crystal in the conditions when 
simultaneously SHG, uj+uj — ► 2uj, and SFM, u+2u; — ► 3ui 
take place. Equations of motion for the slowly varying 
complex amplitudes Ai (I = 1, 2, 3) of the waves arc |l^, 
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where g{z) is a function that equals to +1 (or -1) in 
a single positive (negative) polarization domain of the 
ferroelectric crystal. In this work for the sake of sim- 
plicity, we consider only a periodic alternative domain 
superlattice with a spatial period A. However, g{z) can 
be a quasiperiodic function in the case of nonlinear qua- 
sicrystals ^ ^| . Note that we consider a typical situation 
A < A, where A is a wavelength [|[ ||, @. 

The coupling constants between waves, (3 2 and (3 3 , are 
defined as 

P2,3 = , 

cn 

where d e ff = 2ix~x^ and n is a refractive index. The 
refractive indexes for the different waves are different 
due to a light dispersion. However, it may be shown that 
An/n ~ A/ A -C 1 in the conditions of QPM, therefore in 
what follows we will take (3 2 = (3 3 = f3. Finally, the phase 
mismatches involved in Eqs. ([l]) are Ak 2 = k 2 — 2k\ and 
Ak 3 = k 3 — k 2 — k\. Let both these mismatches be 
compensated by a reciprocal lattice vector of NPC, that 



Ak 2 = 27rmi/A, Afc 3 = 2nm 2 /k, 



where m,- = ±1,±3, ±5, . 



The methods of achieve- 



ment of QPM for several parametric processes in a sin- 
gle NPC have been recently discussed in Refs. || ||, |uj 
(theory) and Refs. [Q, |l3|, O] (experiment). 

The dynamical system (pi) together with initial condi- 
tions, which in our case are the values of complex ampli- 
tudes at the boundary of NPC, Aj(z = 0), completely 
determine the nonlinear spatial evolution of waves. Be- 
fore specification of these initial conditions, we can fur- 
ther simplify the equations of motion. First, we in- 
troduce new scaled amplitudes ai = Ai/(\ // IAq), where 
1 = 1,2,3 and A = max (|Ai(0)|, |A 2 (0)|, \A 3 (0)\). Sec- 
ond, we make the Fourier expansion of the function g(z) 
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ai = -a 2 al - £a 3 a 2 , 
a 2 = 0.5al - £,a 3 a{, 
d 3 = £aia 2 , 



(3) 



where £ = \^m 2 /m 3 (nij are the quasi-phase matching 
orders, see Eq. (|); we assume that m 3 > m 2 ). The 
overdots in Eqs. (|) mean the derivatives in respect to 
z/l n i, and a characteristic nonlinear length l n i is defined 
as 
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In the derivation of motion equations (^) we removed all 
fast varying terms performing averaging over 2n/A. It 
can be shown that such a procedure is correct if /„; A 

The equations (|3j) can be presented in the canonical 
form with the Hamiltonian function 
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Additionally to the energy of wave interaction E = H, 
(Eq. (|^)), the dynamical system (||) has the integral of 
motion 



01 



2|a 2 | 2 + 3|a 3 | 



1 



(6) 



corresponding to the conservation of energy of nonin- 
teracting waves. In a general case the system (||) does 
not have other global integrals of motion, thus it is non- 
integrable and should demonstrate chaotic dynamics for 
many initial conditions a/(0) [jl7| However, for some 
values of £ and some specific initial conditions, an addi- 
tional local integral of motion can arise. Let us list these 
cases because they include the physically important sit- 
uations. 

First, if one of the parametric processes, either SHG or 
SFM, is dominant < 1 or( > 1), then an additional 
integral of motion arises, which is of the Manley-Rowe 
type fl2| . Second, nonlinear dynamics is strongly depen- 
dent on the initial values of two "resonant phases" ^2(0) 
and ^3(0)1 where 



ip 2 = 20i - 2 , Vs = 0i + 2 - 
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and 9j (J 



1,2,3) are the lightwave phases, i.e. aj = 
\dj I exp (-i9j). We found that for ip 2 {0) = ip 3 (Q) = dy- 
namics is always regular. Moreover, using approaches of 
Refs. H |l9), it is possible to show that an additional local 
motion integral exists in this case pp| . In particular, the 
problem of THG (a x (0) = 1, a 2 (0) = a 3 (0) = 0) belongs 
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to this class of initial conditions. Therefore, the spatial 
dynamics of lightwaves at THG is regular (c/Ref. M], 
where an analytic solution has been found). 

We have performed an intensive search of chaotic tra- 
jectories solving the equations of motion (||) numerically 
for two characteristic values of control parameter £ that 
correspond to the experimental situations described in 
works jjj and fl3fl , correspondingly: 

"Set I": The QPMs of first order for both processes, 
mi = 771 3 = 1, £ = \/3 ~ 1.73. 

"Set II": The QPMs of ninth and 33d orders, m x = 9, 
TO 3 = 33, £ = (3\/3)/ll w 0.472. 

We consider several types of initial conditions, which 
cover practically all physically interesting cases (note 
that all these initial conditions satisfy to the restriction 
arising from the integral of motion (§)): 

"Problem 1": a x (0) = a, a 2 (0) 
[1 -a 2 ] 1/2 2- 1 /2 cxp (_^) ; 03 (o) = o, where the 
real parameters <f> and a vary in the ranges —tt < cf> < ir 
and < a < 1, correspondingly. Obviously here 
1-02(0)1 = 1^(0)1 = |0|. 

"Problem 2": oi(0) = [l - 3a 2 ] 1/2 3" 1 / 2 exp (— i6»x), 

a 2 (0) = [l-Sa^V^expH^), a 3 (0) = 
aexp(— 7^3), where —tt < 8j < tt (J = 1,2,3) and 
< a < 3- 1 / 2 « 0.57735. 

"Problem 3": ai(0) = aexp (-£0i), a 2 (0) = 0, a 3 (0) = 
[1 -a 2 ] 1/2 3- 1 / 2 exp(-i6»3), -tt < % < tt (j = 1,3) and 
< a < 1. 

We start our analysis with "Problem 1". This set of 
initial conditions describes, in particular, the THG at 
a = 1 (0 = 0) and the parametric amplification with 
a low- frequency pump at a <C 1 ]l2| ]. In order to in- 
crease an efficiency of energy transformation from a ba- 
sic wave of frequency w to a wave of frequency 3uj, it 
was suggested recently to mix some nonzero signal at 
the frequency 2uj with a basic beam [Q. Such kind of 
initial conditions correspond to a — > 1 (but a^l) with 
different values of phase <fi. 

To distinguish regular and chaotic dynamics we com- 
pute the maximal Lyapunov exponent X ma x for the dif- 
ferent values of initial lightwave's amplitudes, a, and 
phases, 4>. For chaos X max > 0, in contrast X max = 
for a regular motion [jl8|. The dependence of X ma x on 
a for the first order QPMs (set I) is depicted in Fig. |l|. 
For (j> — the initial values of resonant phases, ^2(0) 
and 7/^3(0), are zero corresponding to the integrable limit 
with X m ax — independently on the value of a (not 
shown in Fig. |l|). However, even a small deviation from 
the integrable limit, |^ 2 (0)| = |^ 3 (0)| = \<j>\ = 0.01, re- 
sults in a chaotic motion for a quite wide range of initial 
conditions (dashed line) . Further increase in the value of 
\(j)\ makes chaos more strong (dotted line, \4>\ = 0.1); the 
most strong chaos arises for \(j>\ = 7r/2 (solid line) corre- 
sponding to the initial values of resonant phases |t/'2,3(0)| 
most distant from the integrable limit. 

The motion is always regular for the standard THG 
(a = 1), as well as for some range of a in the vicinity 
of a = 1 (see the right hand side of Fig. nl). A regular 




Figure 1: Dependence of the value of maximal Lyapunov 
exponent on the amplitude of first wave at the boundary 
of optical supcrlatticc a and for the different phases: (/) — 
-it/2 (solid line), = -0.1 (dotted line) and 4> = -0.01 
(dashed line). The first order QPMs (problem 1, set I). 



spatial evolution of lightwaves for a = 1 is shown in the 
upper subplot of Fig. |2j. However, for \<f>\ = ir/2 a strong 
chaos exists already for a ~ 0.95, i.e. for ai(0) = 0.95, 
a 2 (0) w 0.227, a 3 (0) = 0, see lower subplot of Fig. |. 
Thus, the possibility of transition to chaos must be taken 
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Figure 2: Regular (upper) and chaotic (lower) spatial 
evolutions of scaled intensities of lightwaves at the first 
order QPMs. For the upper subplot a — 1 and = 0, 
while for the lower subplot a = 0.95 and <f> = tt/2. 

into account in the application of an additional pump of 
frequency 2ui in order to increase an efficiency of THG 
II- 

We consider now the situation corresponding to the 
left hand side of Fig. [| with a <C 1- This is the para- 
metric amplification with a low-frequency pump Jl2[ | . In 
this case our analysis demonstrates that the evolution 
of waves is weakly chaotic for | ^2,3 (0) j distant from the 
integrable limit. In this regime, the Lyapunov exponent 
has some very small but yet positive value, therefore it 
is very difficult to distinguish a weak chaos from a regu- 
lar motion. In practical terms, it means that one needs 
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to have a very long sample to see the differences be- 
tween regular and weakly chaotic spatial evolutions of 
light waves. 

Now we turn to the consideration of nonlinear dynam- 
ics using the second set of QPM parameters but same set 
of initial conditions ("Set II, Problem 1"). Main results 
on the transition to chaos are depicted in Fig. 0. Again, 
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Figure 3: Same as in Fig. [l] but for the high order QPMs 
(problem 1, set II): (j) — —ir/2 (solid line), 4> = —0.1 
(dashed line) and 4> — —0.01 (dotted line). 

as in Fig. |], 1^2(0)1 = 1^3(0)1 = \<P\ = results in a reg- 
ular motion, while a motion is chaotic for many initial 
conditions if \<f>\ > 0. However, the absolute values of 
Lyapunov exponent are small: Really, maxA mal ~ 0.1 
in Fig. [j], but m&x\ max ~ 0.01 in Fig. ||. Therefore, we 
conclude that multiple interaction of waves employing 
high order QPMs is more stable against a transition to 
chaos in comparison with the case of first order QPMs. 

We consider now a nonlinear dynamics in the case 
when some portion of energy is presented at z = in 
each of interacting waves ( "Problem 2" ) . We present our 
findings in Fig. ||. A strong chaos arises as soon as one 
of the resonant phases becomes different from the inte- 
grate limit 1^2,3(0)1 = (|V>2(0)| = 7T and |^(0)| - tt/2 
for a solid line, |^2(0)| = tt/2 and | -03 (0) | = for a 
dashed line). We should note that for the parameters 
corresponding to a solid curve in Fig. ^ strong chaos ex- 
ists for almost all values of initial wave amplitudes a. 
Chaos is sufficiently weaker for the high order QPMs 
in comparison with the case of first order QPMs: cf a 
dashed line with a dashed and dotted line that corre- 
spond to the same values of phases 9j but to the different 
sets of QPM-parameters. 

Finally, we analyze the set of initial conditions termed 
as "Problem 3" . In particular, it includes the down con- 
version jj^, [l2| or, in other words the fractional conversion 
lu — ► (2/3)ui |I9| ]) in the case a -C 1. For this set of ini- 
tial conditions we did not find visible regions of chaotic 
dynamics. 

In order to reliably distinguish regular and chaotic 
spatial evolutions of lightwaves in conditions of an ex- 
periment, one needs to have many enough characteristic 
nonlinear lengths, l n i, on the total length of the crystal 
L'- L/l n i > 10 j|, [|. Importantly, it appears possible 
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Figure 4: Dependence of the value of maximal Lyapunov 
exponent on the amplitude of third wave at the boundary 
of optical superlattice, a, and for the different phases and 
the QPM orders (problem 2, sets I and II): 6 X = tt/2, 
#2 = 0, 3 = 7T, first order QPMs (solid line); 61=62 = 
9 3 = -tt/2, first order QPMs (dashed line); 6 X = tt/2, 
9 2 = 0, 63 = 7r, high order QPMs (dashed and dotted 
line). 

to meet this condition in the typical NPCs. Really, for 
periodically poled lithium niobate with a period A = 30 
lira, a crystal length L ~ 1 cm, a nonlinear coefficient 
d 33 = 34 pm/V f|, [l3) and a light intensity A% = 0.76 
GW/cm 2 (A = 1.064 /Ltm) [||, we have L/l ni ~ 100. 
Moreover, chaos should be easier observable in the GaAs 
optical superlattice with > 90 pm/V pi[ . 

In summary, we have shown that a simultaneous mul- 
tiwavelength generation in typical nonlinear photonic 
crystals is often chaotic. This fact must be taken into an 
account in realization of compact laser multicolor sources 
for printers, scanners, and color displays based on the 
quasi-phase-matched harmonics generation. 

We should distinguish our results from a recent paper 
ipB) , where nonlinear spatial field dynamics and chaos 
have been studied in a quadratic media with a periodic 
Bragg grating. 

We thank Andreas Buchleitner for discussions and 
support, Martin Fejer for useful comments, and Pekka 
Pietilalncn for critical reading of the manuscript. KNA 
is also grateful to co-authors of his works j|, ^ for a cre- 
ative cooperation. This research was partially supported 
by the Academy of Finland (grant 100487) and the MPI 
PKS. 

References 

[1] S. A. Akhmanov and R. V. Khohlov, Problems of 
Nonlinear Optics, (Gordon & Breach, New York, 
1973) [Original Russian edition: Moscow, VINITI, 
1964]. 

[2] S. A. Akhmanov, and V. G. Dmitriev, and V. 
P. Modenov, Radiotechnika i Elektronika 9, 814 
(1964). 



4 



[3] M. V. Komissarova and A. P. Sukhorukov, Kvant. 
Elcktr. (Moscow) 20, 1025 (1993) [Sov. J Quant. 
Electr. 23, 893 (1993)]. 

[4] K. N. Alekseev, G. P. Berman, A. V. Butenko et 
al. Kvant. Elektr. (Moscow) 17, 42 (1990) [Sov. J. 
Quant. Electron. 20, 359 (1990)]; J. Mod. Opt. 37, 
41 (1990). 

[5] N. V. Alekseeva, K. N. Alekseev, V. A. Balueva et 
al, Opt. Quant. Electr. 23, 603 (1991). 

[6] A. L. Aleksandrovski, A. S. Chirkin and V. V. 
Volkov, J. Russ. Laser Res. 18, 101 (1997). 

[7] O. Pfister, J. S. Wells, L. Hollberg et al, Opt. Lett. 
22, 1211 (1997). 

[8] X. Liu, Z. Wang, J. Wu, and N. Ming, Phys. Rev. 
A 58, 4956 (1998); K. Fradkin-Kashi and A. Ady, 
IEEE J. Quant. Electr. 35, 1649 (1999); S. Saltiel 
and Yu. S. Kivshar, Opt. Lett. 25, 1204 (2000). 

[9] For recent review see, Y. Y. Zhu and N. B. Ming, 
Opt. Quant. Electr. 31, 1093 (1999). 

[10] V. Bcrgcr, Phys. Rev. Lett. 81, 4136 (1998); N. G. 
R. Broderick, G. W. Ross, H. L. Offerhaus et al, 
Phys. Rev. Lett. 84, 4345 (2000). 

[11] J. A. Armstrong, N. Bloembergen, J. Ducuing, and 
P. S. Pershan, Phys. Rev. 127, 1918 (1962). 

[12] A. S. Chirkin, V. V. Volkov, G. D. Laptev, and 
E. Yu. Morozov, Kvant. Elektr. (Moscow) 30, 847 
(2000) [Quant. Electr. 30, 847 (2000)] and refer- 
ences cited therein. 

[13] V. V. Volkov, G. D. Laptev, E. Y.Morozov et al, 
Kvant. Elektr. (Moscow) 25, 1046 (1998) [Quant. 
Electr. 28, 1020 (1998)]. 

[14] M. L. Sundhcimer, A. Villeneuve, G. I. Stegeman, 
and J. D. Bicrlcin, Electron. Lett. 30, 1400 (1994); 
P. Baldi, C. G. Trcvino-Palacios, G. I. Stegeman et 
al, Electron. Lett. 31, 1350 (1995); S. N. Zhu, Y. 
Y. Zhu, and N. B. Ming, Science 278, 843 (1997); 
X. Mu and Y. J. Ding, Opt. Lett. 26, 623 (2000); 
G. Z. Luo, S. N. Zhu, J. L. He et al, Appl. Phys. 
Lett. 78, 3006 (2001). 

[15] Yu. S. Kivshar, T. J. Alexander, S. Saltiel, Opt. 
Lett. 24, 759 (1999); I. Towers, A. V. Buryak, R. 
A. Sammut, and B. A. Malomcd, J. Opt. Soc. Amer. 
B 17, 2018 (2000). 

[16] A. S. Chirkin and D. B. Yusupov, Kvant. Elektr. 
(Moscow) 9, 1625 (1982) [Sov. J. Quant. Electron. 
12, 1041 (1982)]. 

[17] J. Ford and G. H. Lansford, Phys. Rev. A 1, 59 
(1970). 



[18] A. J. Lichtcnberg and M. A. Lieberman, "Regular 
and Chaotic Dynamics" , (Springer, Berlin, 2nd edit 
1991) [Russ. trasl. of 1st edit: Mir, Moscow, 1984]. 

[19] V. V. Konotop and V. Kuzmiak, J. Opt. Soc. Amer. 
B 17, 1874 (2000). 

[20] K. N. Alekseev, 2001, (unpublished). 

[21] C. Zhang, S. X. Yang, Y. Q. Qin et al, Opt. Lett. 
25, 436 (2000). 

[22] O. A. Egorov and A. P. Sukhorukov, Izvestya 
AN Fiz. 62, 2345 (1998) [Bulletin of the Russian 
Academy of Sciences. Physics 62, 1884 (1998)]. 

[23] P. Vidakovic, D. J. Lovcring, J. A. Levenson et al, 
Opt. Lett. 22, 277 (1997). 

[24] L. A. Eyres, P. J. Tourreau, T. J. Pinguet et al, 
Appl. Phys. Lett. 79, 904 (2001). 

[25] A. V. Buryak, I. Towers, and S. Trillo, Phys. Lett. 
A 267, 319 (2000). 



5 



